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The dilaton-gravity sector of a linear in the scalar curvature, scale invariant Two Measures Field 
Theory (TMT), is explored in detail in the context of closed FRW cosmology and shown to allow 
stable emerging universe solutions. The model possesses scale invariance which is spontaneously 
broken due to the intrinsic features of the TMT dynamics. We study the transition from the 
emerging phase to inflation, and then to a zero cosmological constant phase. We also study the 
spectrum of density perturbations and the constraints that impose on the parameters of the theory. 

PACS numbers: 98.80.Cq, 04.20.Cv, 95.36.+X 

As a way to address the cosmological constant (CC) problem [H-[3, the accelerated expansion of the late time 
universef^, the cosmic coincidence ^ (see also reviews on dark energyp-Q, dark matter and references therein), 
many models have been proposed with the aim to find answer to these puzzles, for example: the quintessence [lOj , 



coupled quintessence fc-essencejl2l|.|l^ 

One can add to the list of puzzles the problem of initial singularity |lj|,[l5|, including the singularity theorems for 
scalar field-driven inflationary cosmology pjj, resolution of which is perhaps a crucial criteria for the true theory. The 
avoidance of the initial singularity will be the central question that we will address in this paper, exploring the idea of 
the "emerging universe", where the universe has a non singular origin, such that the Einstein Universe. Although the 
original proposal for the emerging Universe [l7j suffered an instability, several proposals to formulate a stable model 
have been given [l^, in particular one is obtained by invoking Jordan Brans Dicke modelstl9] 

In this paper we explore a model including gravity and a single scalar field (j) in the framework of the so called Two 
Measures Field Theory (TMT)^3|-[1E|- In TMT, many cosmological issues can be addressed: zero vacuum energy 
is obtained without fine tuning, the fifth force problem is resolved and the Einstein's GR is restored when the local 
fermion matter energy density (i.e in the space-time regions occupied by matter) is much larger than the vacuum 
energy density. In this paper, we will address the existence and stability of the emerging universe in TMT. In a 
previous work, [26| in which a square curvature term was almost of the TMT structure (except for a contribution 
that gave rise to a cosmological term), an emergent universe was described. 

TMT is a generally coordinate invariant theory, where the action has to be of the form[2l|-(25| 

j Li<^(i^x + j L2y/^d'^X, (1) 

including two Lagrangians Li and L2 and two measures of integration ^/—g and $. One is the usual measure of 
integration y/—g in the 4-dimensional space-time manifold equipped with the metric g^^. The other is the new 
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measure of integration $ in the same 4-dimensional space-time manifold. The measure $ being a scalar density and 
a total derivative (see Ref. jl^l) may be defined by means of four scalar fields tpa (a — 1, 2, 3, 4), 

^ = e^"'°''^eabcddf,ipadu'fbda'PcdfsiPd- (2) 

It is assumed that the Lagrangian densities Li and L2 are functions of all matter fields, the dilaton field, the metric, 
the connection but not of the "measure fields" (ipa )■ In such a case, i.e. when the measure fields enter in the theory 
only via the measure <1>, the action ([T|) possesses an infinite dimensional symmetry. In the case given by Eq.(l2]) these 
symmetry transformations have the form t^a — ^ </3a + /a(ii), where /q(Li) are arbitrary functions of Li (see details 
in Ref.[21|). In this paper we will insist in keeping this symmetry and therefore the TMT structure. 

We assume here that all fields, including also the metric, connection and the measure fields are independent dynam- 
ical variables. All the relations among them are results of the equations of motion. In particular, the independence of 
the metric and the connection means that we proceed in the first order formalism and the relation between connection 
and metric is not necessarily according to Riemannian geometry. 

Varying the measure fields ipa, we get B^d^Li — where — e^'"^'^ Eabcddufbda'fcdpipd- Since Det{B^) = ^jr^^ 
it follows that if $ 7^ 0, 

Li = sM* = const, (3) 

where s = ±1 and M is a constant of integration with the dimension of mass. 

We proceed now to discuss the question of scale invariance in the context of TMT. A dilaton field (p allows to realize 
a spontaneously broken global scale invariance (22j. We postulate that the theory is invariant under the global scale 
transformations: 



9tj.v ^ e^g,,u, r^^^^r^^, ipa~^\ab^b where det(Aab) = e^^ (t)~>(j)-^e. (4) 



We choose an action which, except for the modification of the general structure caused by the basic assumptions of 
TMT, does not contain any exotic terms and fields as like as in the conventional formulation of the minimally coupled 
scalar-gravity system. Keeping the general structure ([1]), it is convenient to represent the underlying action of our 
model in the following form [2J|: 



S= J d\e'^^/'"^ __Li?(r,g)(cI> + hg^g) + ($ + b^^g)^g^''^,^cj,^, - e'^^"'^ + ^gV^) 



(5) 



where bg represents the coupling constant of the curvature scalar and of the scalar kinetic term to ^ —g, respectively. 

We use K = ^tt/M"^ where Mp is the four-dimensional Planck mass. In the equations of motion following from this 
action, we change the metric to the new one 

5^, = e«^/^p(C + 6,)g^„ (6) 

where C, = . The conformal metric represents the " Einstein frame" , since the connection becomes Riemannian. 
Notice that g^^, is invariant under the scale transformations (|4]). After the change of variables to the Einstein frame 
the gravitational equations take the standard GR form 

G^,{g^p) = KT;fJ, (7) 
where G^j,u{gai3) is the Einstein tensor. The energy-momentum tensor, T^l^ , becomes 

= 1^ - \g,.r^<^,c.ct>j^ - g^- sCC + M ^"'''^'"'^'^ + 5M^'^e//(0; C, Af), (8) 

where the function Vef f{<p]C,,M) is defined as following: 

K//(</';C,M) = ^i ^^^^^^2 ^ (9) 



where Vi and V2 are two arbitrary constants. In order to have 14// > the two constant, Vi and V2, should satisfy 
the inequality bgVi > V2 for — > 00. 
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The scalar field ( is determined by the consistency of ([7]) with ([3]), which lead to the constraint 

{bg ~ C) [sM^e-^"^/*^^ + Vij -2V2-6- 6g(C + bg)Z = 0, (10) 

where Z = Ig'^f^ (jj^act),^ and 5 = 

The effective energy-momentum tensor ([8]) can be represented in a form of a perfect fluid T^l^ — {p+p)u^u^ -~P9tii^, 
where = {2zy-/2 with the following energy and pressure densities resulting from Eqs.® and ^ after inserting the 
solution C = C(0, Z; M) of Eq.([TU)) 

(sAf4e-2"'^/^^P + 14)2 _ 2,56„(sM4e-2«'^A'^p + \/i)Z - S-S^b^Z^ 

4feUe-t/M.^v..)_v.,] (") 

and 

^ A[bg{sM'^e-^°"l'/^^^ +Vi)-V2] ^ ' 

Notice that if s and V\ have different signs one might obtains a state with zero energy density. 

We now want to consider the detailed analysis of The Emergent universe solutions and in the next section their 
stability in the TMT scale invariant theory. We start considering the Friedmann- Robertson- Walker closed cosmological 
solutions of the form 

ds^ ^ - a{tf + r'^{Ad'^ + sin^Bd<\?)^ , (f) = (t>{t), (13) 

where a(t) is the scale factor, and the scalar field (f) is a function of the cosmic time t only, due to homogeneously and 
isotropy. We will consider a scenario where the scalar field (f) is moving in the extreme right region — > oo. In this 
case, the expressions for the energy density p and pressure p are given by, 

P = + SB^'' + C, (14) 

and 

p^^^^ + B^^-C, (15) 
respectively. Here, the constants A, B and C are given by, 

A = l- , ^ ' B = 2 -, and C=- -, (16) 

^{bgV^-V2)' ^{bgVi-V2)' ^{bgVi-V2)' ^ ' 

It is interesting to notice that all terms proportional to (jr' behave like "radiation", since p^i = is satisfied. In 

the same way, the terms proportional to cfP' behave like "stiff' matter, since p^^ = p^2, and finally, the C constant 
term behaves like a "cosmological constant". The emerging universe can turn into inflation only if C > 0. 

The equations that determines the static closed universe a{t) = ao = constant, in which a — 0, d — 0, gives rise to 
a restriction for (jjQ that have to satisfy in order to guarantee that the universe be static. Since, a = is proportional 
to p -f 3p, we must require that p + 3p = 0, which leads to 

3B04 + A<j>l-C = 0. (17) 

This equation leads to two roots, given by 



respectively. 



.,^^±VArTT^-A__ ^^^^ 



6B 
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Defining the variable y = we see that A = 1 — y, BC = — fg, so the condition that the discriminant be 

positive, i.e., that A^ + UBC > gives (1 - y)^ - 2|i > or ^ - 2y + 1 > 0, which is satisfied for y < 2(2-^3) = 0.54 
or y > 2(2 + -s/S) — 7.46. However, for C > the constraint (/)^ > is only satisfied for y < 0.54 and the other region 
y > 7.46 has to be discarded. Stability provides further constraints. In fact, as we will see, the second solution can 
never be stable. 

It is also interesting to see that if the discriminant is positive the first solution yields automatically a positive energy 
density, if we require C > 0. The same requirement is to be adopted if we want the emerging solution to be able to 
turn into an inflationary solution. Since C > 0, we get that {bgVi — V2) > in agrement with the fact that Veff > 
as was mentioned previously, and therefore we get that B < 0. One can see that the condition p > for the flrst 
solution reduces to the inequality w > {1 — \/l — w)/2, where w = —12BC/A^ > (recall that B < 0), and as long 
as w < 1, it is always true that this inequality is satisfied. 

In the following we will study the stability of the static solution. Let us consider the perturbation equations. 
Considering small deviations for from the static emergent solution (^o) and also the perturbation of the scale factor 
a, from Eq. ([Ti)) we obtain that 

dp ^ ^00(50 + l2B(i)ld(j). (19) 
At the same time 6p can be obtained from the perturbation of the Friedmann equation, i.e. 



1 



3 ( ^ + ^ ) = '^P' (20) 



and since we are perturbing a solution which is static, i.e., H = 0, we obtain that 



—Sa = K5p. (21) 



We also have the second order Friedmann equation, 

1 + + 2ad 



2 - -^^P■ (22) 

Applying to this equation the static emergent solution, i.e. po = — Po/S and a = ao, we get 

2 2 

— = -2kpq = -Kpo = JloK/Oo, (23) 
Oq 6 

where we have chosen to express our result in terms of flo, defined by po = [Qq — l)po, which for the emerging solution 
has the value fio = |. Using this in Eq. ipTt . we obtain 

Sp^-^Sa, (24) 
ao 

and equating the values of Sp as given by Eqs. ([T9l) and (f24|) we obtain a linear relation between d(j) and 5a, which is 
given by 

S<j) = DoSa, (25) 

where 

^ . 3^"Po ^2g) 

aoMA + l2Bcj)l) 

We now consider the perturbation of the Eq. (f22|) . In the right hand side of this equation we take that p = (fi — l)p, 
with 

n^2(l-^), (27) 
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where 

U,ff = C + B^\ (28) 
and thus, the perturbation of the Eq. (|24| leads to, 

-'^ + 2—^-K6p^^Kdm-l)p]. (29) 
% ao 



In order to evaluate this, we use Eqs. ([27| and (|28p . and the expressions that relate the variations in a and (j) given 
by Eq.([25l). Defining the "small" variable /3 as (/3 < f) 



a(i) = ao(f + /?), 



(30) 



we obtain 



(31) 



where 



w^ = nopo 



A- 



12 (j)lB 



6 



Po 



Snflo + 2k 



(32) 



Notice that the sum of the last two terms in the expression for W^, that is SkUq + 2k vanish since fto — 



the same reason, we have that 6 



{C+B4>t) _ 



4, which brings us to the simplified expression 



|. For 



6B 



A +12 



blB 



(33) 



For the stability of the static solution, we need that Wq > 0, where 0o is defined either by Eq. (0o = 'Pi) oi" 
— v^i)- Notice that since B < 0, 2A B (pQ < 0, so that for W§ > 0, we need A + 12 0§ S < 0. Wg > implies 
> 1. Multiplying this inequality by the negative quantity A + 12 (PqB < and evaluating separately for 



6B 



A+12<j>l B 



02 ^ 02 g^j^j ^2 ^ 02 ^^^^ ^Yie condition Wq > becomes iy/A"^ + 12BC < Sy/A"^ + 12BC if we use the first 

solution and —Ay/A'^ + 12BC < —8vA^+12BC if we use the second solution. Of course this means that the second 
solution can never be consistent with stability. For the first solution, we still have to verify that A + 12 (fy^B < 0. 



- 

2Sb„C 



Vi 



Introducing the relevant expression for 0o appropriate for the first solution, and using again the variable y 
obtain now that y > 1/2. Putting all together, the existence, the stability previously defined, as well as a reasonable 
cosmological picture after the emerging phase (inflation), are satisfied for y < 1, to avoid negative kinetic terms during 
the slow roll phase of inflation. Therefore, in order to have a picture in which the emerge universe is stable and then 
pass to an inflationary phase is obtained if the following range is satisfled 



0.5 <y < 0.54. 



(34) 



The study of the stability of the static solution and the properties of the different equilibrium points could be done 
in a more systematic way by using a dynamical system approach. In this scheme we rewrite the Friedmann and the 
conservation of energy equations as an autonomous system in terms of the variables H = d/a and x = (fP' . In order 
to do so, we differentiate the Fridmann equation and after using the expressions for p and p, given by Eqs. (ll4l) . (jlSp 
and (1221). we obtain: 



3 



C -iBx^ - Al 



(35) 



and 



-6H 



Ax + ABx"^ 
A+12Bx ' 



(36) 
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where A = 1 — y, as before. The equations (j35|) and p6p are a two-dimensional autonomous system on the variables 
H and x. 

In order to study the stability of the static solutions we look for critical points of the system (|35|) and (|36| . These 
points are 



^ff = 0, X 



-A + y/A^ + 12BC y 



6B 



{ 



iJ = 0, x = 



-A-VA^ + ^2BC^ 



6B 



{h-~J^,^ = 0};{h=J^,x^0};{ 



H 



{ 



{A^ + 16BC)k 
B 



4^/3 



H = 



(A^ + 16BC)t, 



4^/3 



-,x 



4B/' 

= --}■ 
ABi 



(37) 
(38) 

(39) 
(40) 



The critical points have different properties depending on the values of the parameters of the model {B and C). At 
this moment we are not going to give an exhaustive description of these properties for all the critical points, instead, 
we are going to focus on the particular critical points which are related with static universe. From the definition of 
the variables H and x we can note that only the first two critical points Eqs. (157)) and ([55]) correspond to a static 
universe. In order to study the nature of these two critical points we linearize the equations psp and p6p near these 
critical points. From the study of the eigenvalues of the system we found that the first critical point, Eq. ([57)1 . could 
be a center or a saddle point, depending on the values of the parameters of the model. On the other hand, the second 
critical point Eq. ([38]) is a saddle. 

Stable static solutions correspond to a center. This imposes the following conditions for the parameters {B and C) 
in order that the critical point, Eq. (j37|) . becomes a center. 



B<0, <C<-i/^ . (41) 

' 64B y B^ AB ^ ' 

These conditions also ensure that x > and the positivity of the energy density. If we consider the definition of 
the parameters A, B and C given in Eg. dTB]) we can note that Eq. (HIT) is in agremcnt with the stability conditions 
that were found previously. 

In Fig. [Hit is shown a phase portrait near the center critical point {H — 0,x — 0.0443) for three numerical solution 
to Eqs. (1551) and ([55)) . Also, in this figure we have included the Direction Field of the system in order to have a 
picture of how a general solution look like. In this figure we have used the values B — —1 and C = 0.016. 



♦ fit** 

f f f f f f f * , 
t f t f f t f t f 



^ ol)42t' 'kcM-to'^ o^d435^ 'll.olim '-elof 15 r o!o45( 




t t t t 1 t t * 

» t t ♦ 



FIG. 1: Plot showing the Direction Field near the center critical point and three numerical solution. Here we have used unit 
where n — 1. 
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Once a transition to a slow roll phase inflationary phase takes place, we then have to see if the resulting inflationary 
phase can provide enough e-foldings for the solution of the Big-Bang standard problems. 

We consider the relevant equations in the slow roll regime, i.e. for 0^/2 ^ V{(t)) and when the scalar field is 
large, but finite. Dropping higher powers of (j) in the contributions for the kinetic energy, we obtain. 



1 



with 



and 



7(0) 



eff^ 



4[6<,(M4e-2«'A/^^p +Vi)-V2y 



^bgVl - V2] 



(bgVl-Va) 



(42) 
(43) 

(44) 
(45) 



Here we have taken s — \, hg > Q and Vi < 0, then one obtains without fine tuning a vacuum state with zero energy 
density and thus Veff{4> = 0, </>) > 0. 

In the slow roll approximation, we can drop the second derivative term of and the second power of (j) in the 
equation for and obtain. 



(46) 



where V^jj: — -^j^- The relevant expression for Vefj will be that given by Ea. (|^5)) . i.e., where all higher derivatives 
are ignored, consistent with the slow roll approximation. 

We now display the relevant expressions for the asymptotic value of 0, these are 



(47) 



where 



and 



Vo 



^{bgVi-V2) 



V 



Vl M\bgVl ~ 2V2) 

^{bgVi-V2Y 



2{bgVi - V2) 



-2a4i/Mj, 



SbgV2M^ 

2{bgVl ~ V2f 



70 + 71 e 



(48) 



Note that F < since 6g > 0, \4 < and V2 < 0. 

At the end of inflation, where, (/) — 4>end^ the parameter e, defined by e = —jp-, takes an approximated value equal 
to one (analogous to a ~ ) . The condition under which inflation takes place can be summarized with the parameter 
£ satisfying the inequality £ < 1 (or d > 0). Taking the derivative with respect to the cosmic time of the Hubble 
parameter and from Eq. pSl) . we obtain that the condition e ~ 1 gives 



- = ^(k:///k//)^-i, 

working to leading order, setting 7 « 70, K// ~ Vb and V^^j: ~ — {2a/ Mp) Vexp{—2a(l)/Mp), we obtain 

V2a \ V\ ' 

We now consider 0* and the requirement that this precedes (j)end by N e-foldings. 



N : 



Hdt ; 



—I 



Kff 



(49) 



(50) 



(51) 
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In the following, the subscripts * and end are used to denote to the epoch when the cosmological scale exit the horizon 
and the end of inflation, respectively. 

Solving in terms of Veff using Eq. (|46p . working to leading order, setting 7 = 70 and integrating, we obtain that 
the relation between and N becomes 




V2- 



(52) 



We finally calculate the power of the primordial scalar perturbations. The power spectrum of the curvature per- 
turbation in the slow-roll approximation for a not-canonically kinetic term becomes Ref.[2^(see also Refs.[29|) 



£j2 

Ps = ki , (53) 

Cs e 

where it was defined the "speed of sound", c^, as = p z+izp zz ' ^^^^ P{Z^'I>) function of the scalar field (j> and 
the kinetic term, Z = {l/2)g^'''^d^(j)dv4>, and fci = (87rMp)~^. Here P z denote the derivative with respect Z. In our 
case P{Z, (j)) = 7(</>) Z -Veff, with Z = <p^/2. Thus, from Eq.® we get 

Ps = ki ^. (54) 

7(0)02 

The scalar spectral index is defined by 

".-l = ^^ = -2£-^-^, (55) 

where rj — ^jj ^^d f — -^^^ . 

On the other hand, the generation of tensor perturbations during inflation would produce gravitational waves. The 
amplitude of the tensor perturbations was evaluated in Ref . (28j . In our case 

2 f 2ZP^z-P \ 
^^ = ^ [ —i ) ' (5*^) 



37r2 V Mp y 

where the tensor spectral index, ut, becomes nx — '^J"^^^^ — ~2e, and they satisfy a generalized consistency relation 
given by r = = —8 nx- 

Therefore, the scalar field (to leading order) that should appear in Eq. (|54| should be -^700 and thus we have 




Vo fcl 

12 



r- AaN 

V2- 



2 



(57) 



This quantity should be evaluated at = 0* given by Ea. dS^ . 

In Figl2]we show the dependence of the tensor-scalar ratio r on the spectral index n^. From left to right a = 0.01 
(solid line), a — 0.05 (dash line) and a = 0.1 (dots line), respectively. The dots represent the number of c- folds 
for the value N = 60. From Ref. [13], two-dimensional marginalized constraints (68% and 95% confidence levels) on 
inflationary parameters r and ris, the spectral index of fluctuations, defined at kg = 0.002 Mpc~^. The seven-year 
WMAP data places stronger limits on r. In order to write down values that relate Us and r, we used Eas. ([55)) . ([55)) 
and (|57|) . Also we have used the values y = 0.52, M = 1 and k = 1, respectively. We noted that the parameter a, 
which lies in the range 1 > a > 0, the model is well supported by the data as could be seen from Figl2l 

The dilaton dependence of the effective Lagrangian appears only as a result of the spontaneous breakdown of the 
scale invariance. If no fine tuning is made, the energy density p{(j), Z) and the pressure p(0, Z) depend quadratically 
upon the kinetic term Z. Hence TMT represents an explicit example of the effective k- essence resulting from first 
principles without any exotic term in the underlying action intended for obtaining this result. These non linearities 
in p(0, Z) and p(0, Z) play a crucial role in existence and stability of the emerging universe solutions. In this paper 
we have been successful in describing an emergent universe in a TMT sort of theory. 
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0.4 



0.3 



0.2 



0.1 



0,0 

0.92 0.94 0.9() 0.9i> 1.00 1.02 



FIG. 2: The plot shows r versus Us for three values of a. For a — 0.02 solid line, a — 0.05 dash line and a = 0.1 dots line, 
respectively. Here, we have fixed the values y = 0.52, M = 1 and k = 1, respectively. The dots represent the number of e-folds 
for the value N = 60. The seven-year WMAP data places stronger limits on the tensor-scalar ratio [30^ . 
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